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Abstract 

Density functional theory is utilized to articulate a particular generic deconstruction of the electrode/electro-catalyst assembly for 
the cathode process during water splitting. A computational model was designed to determine how alloying elements control the 
fraction of H2 released during zirconium oxidation by water relative to the amount of hydrogen picked up by the corroding alloy. 
This model is utilized to determine the efficiencies of transition metals decorated with hydroxide interfaces in facilitating the 
electro-catalytic hydrogen evolution reaction. A computational strategy is developed to select an electro-catalyst for hydrogen 
evolution (HE), where the choice of a transition metal catalyst is guided by the confining environment. The latter may be recast into 
a nominal pressure experienced by the evolving H2 molecule. We arrived at a novel perspective on the uniqueness of oxide 
supported atomic Pt as a HE catalyst under ambient conditions. 



Introduction 

Molecular hydrogen produced by water splitting constitutes the 
archetypical energy carrier in chemistry and is a main target 
process for the future harvesting of solar energy. Today, water 
splitting represents large economical values, i.e., it comprises a 
significant fraction of the total industrial electric energy 
consumption in the USA [1]. Decisive factors jointly deter- 
mining the efficiency of the electrochemical process are the 
reactions at the oxidizing anode as well as at the hydrogen 



evolving cathode. In two inspiring experimental studies [2,3], 
Subbaraman et al. reported enhanced hydrogen evolution 
activity in water splitting by tailoring TM(OH)2-Pt electro-cata- 
lyst/electrode assemblies, where TM represents Mn 2+ , Fe 2+ , 
Co 2+ and Ni 2+ . The role of these hydroxides was to catalyze 
water dissociation. In this context, the objective of the present 
study is to contribute a novel descriptor for the electro-catalytic 
hydrogen evolution reaction (HER). It offers a complementary 
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perspective on a recent study addressing the oxidation of zirco- 
nium alloys by water [4,5]. The overall reaction 



l A Zr + H 2 0 ->■ Vi Zr0 2 + H 2 



(1) 



is taken to occur by water utilizing hydrolysis to penetrate 
the oxide scale along hydroxylated grain boundaries, see 
Figure la, 



Indeed, Equation 6 was found to be decisive for the fraction of 
hydrogen atoms not forming H 2 but being absorbed in the Zr 
alloy according to 



^HO-Zr IV -H-M X J+ aZr 
[Zr™ - O - M X J + Zr a H 2 (s). 



(7) 



[zr IV -0-Zr IV l + H 2 0 

L Jgrain boundary 

[zr IV -0-H HO-Zr IV ] 

L J; 



(2) 



hydroxylated grain boundary 



These hydroxide ions subsequently react with transition metal 
decorated sites (see Figure la) and zirconium metal to produce 
ZrC>2 in conjunction with transient transition metal associated 
hydride-proton (hydroxide) pairs (see Figure lb) to restore the 
ZrC>2 grain boundary according to 



Zr IV -0-H HO-Zr IV J + | Zr IV -0-M X l + 'AZr 



• [zr IV -0-Zr IV 



] + [HO - Zr IV - H - M X ] + '/2Zr0 2 



(3) 



This can be subdivided into an anode process 



ViZx + 



Zr IV -0-H HO-Zr IV 



-> Vi ZrO, + Zr IV - O- Zr™ + 2H + + 2e~ 



(4) 



For completeness, a 1.1 eV/H2 drive to release H2 from the 
confining grain boundary was computed according to 



[H 2 ] 



grain boundary 



"> H 2(g) + [] 



grain boundary ' 



(8) 



Utilization of the hydride-proton recombination channel (see 
Figure lb), the correlation between the computed reaction ener- 
gies for the HER, Equation 6, and the experimental hydrogen 
pick-up fractions (HPUF), i.e., the fraction of the hydrogen, 
which does not undergo hydrogen evolution but are instead 
picked-up by the alloy during zirconium oxidation by water, 
leads to a model as displayed in Figure le. It is noteworthy, that 
the energetics for the chemical reaction step in Equation 6 
offers a measure of the confinement-dependent cathodic over- 
potential for the HER along the reaction channel (Equations 
2-6). The relevance of the reversible hydride-proton recombina- 
tion reaction (Equation 6) has recently been proposed in case of 
a nickel electro-catalyst supported by seven-membered cyclic 
diphosphine ligands containing one pendant amine, with the Ni 
supporting the hydride and the amine providing the proton in 
the hydride-proton recombination reaction [6]. 



where the [Zr IV -0-Zr IV ] oxide grain boundary is recovered, 
and a cathode process 



Zr IV -0-M X 



Zr IV_ 0 _M X - 2 



HO-Zr IV -H-M 



+ 2H + + 2e~ 

+ 2H + 
x" 



(5) 



is employed to decide the oxidation state X. The subsequent 
chemical drive for H2 release into the confining grain boundary 
determines M and recovers the [Zr IV -0-M x ] site (cf. 
Figure lc) 



[HO-Zr^-H-M 5 ^ 



-H Zr IV -0-M X 



] + [H 2 ] 



(6) 



'grain boundary 



Results and Discussion 

In the following we introduce and employ the notion of 
"confinement effect" as a steric Pauli repulsion type interaction 
between H 2 and a hydroxylated interface (see Figure lc) upon 
hydride-proton recombination. First, we employ this notion in 
the context of hydride-proton recombination reactions to 
demonstrate how it decides which oxidation state X of metal ion 
M minimizes the overpotential for the HER, as quantified by 
the reaction Equation 6 (cf. Figure ld,e). Second, it is shown 
how the emerging understanding is naturally extended to 
include electro-catalysts for HER under ambient conditions. 

Impact of confinement on HER during zirco- 
nium oxidation by water 

To investigate the confinement effect on the HER, we consider 
the zirconium oxidation by water (see Figure Id). The differ- 
ence between the two horizontal lines at 2.9 eV and at -0.2 eV 
corresponds to the 3.1 eV/H 2 [9] thermodynamic drive for HER 
in case of Zr oxidation by water under ambient conditions 
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Figure 1 : (a) Representative structure for a model of a hydroxylated inter-grain interface comprising ZrO(OH>2 on Zr02. Here, this interface is deco- 
rated with Fe ions with the oxidation state +2. Oxygen is represented as red, zirconium as light blue, iron as purple, and hydrogen as white, (b) One 
hydride ion and one hydroxide moiety prior to the hydride-proton recombination to form H2 is displayed, reactant in Equation 6. (c) The product in 
Equation 6 is displayed, including M x coordination to the additional oxygen ion replacing the hydride ion and the released grain boundary H2. (d) 
Hydride-proton recombination energies for H2 release into said interface (dashed black line at 1 .1 eV), enthalpy change for H2 release at ambient 
pressure (dashed black line at 0 eV), and corresponding Gibbs energy change (dashed black line at -0.2 eV). TM 2+ blue, TM 3+ red, weighted average 
green, (e) Comparison of theoretical data (green) and experimental HPUF data (black); * from [7] and 0 from [8]. The theoretical data is a weighted 
average between TM 2+ and TM 3+ . The black dashed line is HPUF in pure Zr02 from [7]. The blue dashed line corresponds to HE from Zr 4 " 1 " hydride at 
GB with Na + , Ca 2+ and Sc 3+ spectator. Sc 3+ corresponds to the top line, Na + to the middle line and Ca 2+ to the bottom line. 
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(Figure 2). The line at 0 eV represents the energy of a free H2 
molecule at 0 K. The line at 1.1 eV represents the energy cost at 
0 K for bringing a free H2 molecule into the confinement repre- 
sented by Figure lc. The line at -0.2 eV is owing to the 
increase in entropy when a water molecules H2CKI) is consumed 
(-70.0 Jmor'K" 1 [10]) and a H 2 (g) molecule is formed 
(130.7 JmoP'K" 1 [10]) at 298 K and 100 kPa (compare Equa- 
tion 1), while neglecting changes in entropy in Zr upon oxi- 
dation. 

Thus, a perfect electro-catalyst would exhibit an enthalpy 
change of 0 eV for the HER under ambient conditions. More- 
over, it is inferred that a perfect electro-catalyst, which passes 
the HER into this hydroxylated interface via Equation 6 prior to 
the subsequent H2 release under ambient conditions, must 
display 1.3 eV/H 2 overpotential, i.e., (1.1 - (-0.2)) eV/H 2 . 
Equivalently, in case of the HER into the interface, any residual 
drive towards H2 formation relative to the line at 1.1 eV/H2 
corresponds to a local overpotential for the HER into the 
confining interface. A correlation emerges between a greater 
overpotential and a lower hydrogen pick-up fraction (HPUF, 
see Figure Id and le). Thus the well-known effect of Ni 2+ to 
cause detrimental hydrogen pick-up was explained by its reluc- 
tance to release H2 into the hydroxylated internal inter-grain 
interface [4,5]. From the overall agreement between reaction 
energies for Equation 6 and the experimentally reported 
HPUF's, it was concluded that "anti-catalysts" are preferred in 
order to mitigate the HPUF. In case of zirconium oxidation by 
water, these "anti-catalysts" are ions, which conserve signifi- 
cant fractions of the drive for hydrogen evolution by forming 
highly reactive metastable hydrides. These species are 



contrasted by Co 2+ and Ni 2+ , which catalyze the HER when H2 
is released into the highly constraining interface (see Figures Id 
and le). 

A stability check on the semi-quantitative validity of the model 
was offered by a comparison of the experimental 44% hydrogen 
pick-up fraction of pure zirconium (corresponding to the black 
horizontal dashed line at 1.5 eV/H2) and model calculations for 
the hydride-proton recombination energetics employing the 
inert Na + , Ca 2+ , and Sc 3+ as "dummy" ions in the positions of 
the transition metal ions (see the three blue horizontal dashed 
lines in Figure le). 

On HER at ambient conditions - a consis- 
tency check 

According to the above understanding, which ions constitute 
viable electro-catalysts in the absence of confinement or at 
atmospheric pressure? In as much as the drive for HER 
comprises the relaxation of the resulting oxy-hydroxy ions coor- 
dinating the transition metal ion [4], it is suggested that besides 
being able to form the hydride intermediate, metals with low 
oxidation states and large ionic radii should be considered in 
order to minimize their affinities to the oxide surrounding. This 
characterization clearly points to the noble metal ions as candi- 
dates for electro-catalysts. Additional requirements for any 
successful electro-catalyst include sufficient electron conduc- 
tivity of the oxide matrix supporting the catalyst as well as elec- 
tric contact to the electrode itself. Finally, the "water dissocia- 
tion catalysis" put forward by Subbaraman et al. [2,3] is used to 
infer that hydroxylated interfaces provide natural channels for 
proton transport to the oxy-hydroxide supported electro- 
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Figure 2: The diagram on the left is identical to Figure 1d. The enlarged region exposes the overpotentials for the elements in the Pt and Au groups. 
Note that Pt + associated hydride displays a negative overpotential implying that it is more stable than the H2(g) asymptote (lower dashed line). 
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catalytic site. A schematic representation of this understanding 
of the electrode/electro-catalyst assembly is provided in 
Figure 3. 

Employing the above described hydroxylated inter-grain inter- 
face model as a generic supporting matrix for the eletrocat- 
alytic process, we evaluate the energetics for the hydride-proton 
recombination reaction and arrive at a possible descriptor for 
the HER, which is also applicable under ambient conditions. 
This facilitates a procedure for the screening among candidate 
electro-catalysts. 

Indeed, in Figure 2, a descending staircase-like curve for 
electro-catalysts is arrived at for the reaction energy corres- 
ponding to Equation 6. Starting at the hydroxylated zirconia 
inter-grain confinement, where Co 2+ and Ni 2+ are the obvious 
candidate catalysts, we approach the ambient conditions step by 
step by considering the embedded Ag + , Ni + , Cu + , Pd + , Au + and 
eventually Pt + . The Sabatier principle applies in two ways. 
Firstly, X in M x can be made to satisfy the requirement that the 
reactant in Equation 6 forms spontaneously [4,5]. Secondly, the 



environment confining H2 in Equation 6 can be tuned to 
equalize the stabilities of reactants and products, so that any 
drive to release H2 into a confinement is balanced owing to the 
replacement of the three-centered hydride by an oxy-bridge 
(Equation 6). It is gratifying to find that the oxy-hydroxide 
supported Au + and Pt + sites become preferred only when 
approaching ambient conditions, as oxides of noble metals are 
generally unstable, a property which is often associated with 
their softness. Consequently, the drive to replace the hydride 
ion by an oxygen ion is weak, and hence the H2 release is 
expected to require a loose confinement for these systems. This 
is in contrast to harder ions, which form more stable oxides. 
Interestingly, in case of Pt + , the hydride comes out more stable 
than the H2(g) limit. This implies that the embedded Pt + site 
could constitute an efficient absorber of H2 in the gas phase 
under ambient conditions - a purely chemical property. The 
semi-quantitative nature of the methodology does not allow for 
precise predictions of absolute numbers (see horizontal "error 
bar" in Figure le). However, it may be that the overpotentials 
reported for the Pt-based catalysts are related to the coverage 
dependence of the electrochemical decomposition of the Pt + 
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Figure 3: HER at electro-catalyst/electrode assembly. (A) Coalescence of proton and electrons to form the metal catalyst (MC) associated hydride, 
(B) Hydride-proton recombination to form H 2 at the interface. (C) The step between panel B and panel C comprises the HER following the hydride- 
proton recombination step. 
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associated hydride compound. Detailed properties of the 
embedding materials (e.g., electron conductivity) could cause 
the additional variations of the overpotential observed by 
Subbaraman et al. [2]. 

Interestingly, +2 is not considered a relevant oxidation state in 
case of Pt under ambient conditions for the electro-catalytic 
reaction path involving the hydride-proton recombination reac- 
tion (see Figure 2). This result is due to the strong binding of +2 
to the oxy-hydroxide ligands upon H2 release, violating 
Sabatiers principle. 

In conclusion, the present approach offers a complementary 
computational strategy to rank catalysts for HER from water 
splitting. The complex modelling of heterogeneous HER 
electro-catalysis at the interface between composite catalyst/ 
support and a water based electrolyte is subdivided into a chem- 
ical oxide hydrolysis step (Equation 2), an electro-chemical 
redox step (Equation 3, Equation 4 and Equation 5), followed 
by the chemical hydride-proton recombination step 
(Equation 6). This conceptual deconstruction aims for 
supporting the prediction of novel approaches to improve on 
existing electro-catalyst/electrode assemblies. Thus, the design 
of the aqueous electrolyte/substrate system impacts only the 
hydrolysis step (Equation 2). The oxidation state X of M x is 
decided by Equation 5, while the choice of supported HER cata- 
lyst M x is determined according to Equation 6 by the confine- 
ment effect in conjunction with Sabatier's principle. 

For the HER step, a recently proposed alternative to the 
Volmer-Heyrovsky mechanism was employed [4,5]. Rather 
than electron-proton discharge over an M-H moiety resulting in 
the conversion of 2H into H2, the HER investigated here results 
from a hydride-proton recombination reaction. While the 
protons constitute hydroxides in Equation 6, which is non- 
significant due to their ubiquity in aqueous media, an observa- 
tion of three-center hydride intermediates is the sought-after 
"smoking- gun" evidence for the proposed mechanism. 

Computational details 

The Perdew-Burke-Ernzerhof generalized gradient approxima- 
tion PBE GGA [11] as implemented in the DMOL3 engine 
[12,13] in the Material Studios program package [14] was 
employed in conjunction with a double-^ numerical basis set 
with an extra polarization function on each heavy atom and a 
p-function on each hydrogen atom. Systematic spin polarized 
calculations were performed. A 4 x 4 x 1 k-point set for 
sampling the Brillouin zone was compared to a 2 x 2 x 1 
k-point set, and the latter was found to suffice. In order to 
reduce the computational effort, inert electron shells were 
described effectively by means of the semi-core pseudopoten- 



tials. Zero-point corrected free energies were compared to non- 
corrected reaction energies and the differences were deemed 
negligible. 

The grain boundary model (cf. Figure la) was constructed by 
inserting one unit cell of monoclinic ZrO(OH)2 (5.4 Ax 10 A x 
5.4 A) in between two supercells of monoclinic Zr02 (5.4 A x 
10.8 A x 5.4 A), where the unit cell doubling is in the b-direc- 
tion. The stability of the model has been extensively investi- 
gated, including full geometry optimization, when arriving at 
the foundation of [4]. The choice of the grain boundary model is 
far from unique. Here, it is the success in reproducing the 
experimental volcano shape curve (cf. Figure Id) which renders 
the present investigation meaningful. 

The grain boundary model employed to evaluate the reaction 
energy of Equation 6 was subjected to periodic boundary condi- 
tions, where the studied super-cell contained approximately 
50 atoms. The number of hydrogen atoms, i.e., hydroxides and 
hydride, varied. This was because the oxidation states of the 
transition metal ions were controlled by adding (removing) 
hydrogen atoms to (from) the super-cell. This way, neutral 
super-cells were employed in all cases. When evaluating Equa- 
tion 6, all bond distances and bond angles associated with atoms 
in the super-cell were optimized, while the super-cell dimen- 
sions were kept constant. 
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